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We suggest a governing equation which describes the process of polymer chain translocation 
through a narrow pore and reconciles the seemingly contradictory features of such dynamics: (i) a 
Gaussian probability distribution of the translocated number of polymer segments at time t after 
the process has begun, and (ii) a sub-diffusive increase of the distribution variance A(t) with elapsed 
time, A(t) oc t a . The latter quantity measures the mean-squared number s of polymer segments 
which have passed through the pore, A(t) = ([s(t) — $(t = 0)] 2 ), and is known to grow with an 
anomalous diffusion exponent a < 1. 

Our main assumption - a Gaussian distribution of the translocation velocity v(t) - and some 
important theoretical results, derived recently, are shown to be supported by extensive Brownian 
dynamics simulation which we performed in 3D. We also numerically confirm the predictions made 
in ref . (141] . that the exponent a changes from 0.91 to 0.55, to 0.91, for short, intermediate and long 
time regimes, respectively. 
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^ ■ I. INTRODUCTION 



Polymer translocation has attracted recently a lot of attention as it plays a crucial role in processes such as DNA 
and RNA transport through membrane channels [l| . Ultimately, this research might lead to a longstanding objective 
of nucleotic transport: fast and cheap sequencing of DNA. Meanwhile many interesting experimental and theoretical 
findings have been reported |2| . A lot of important observations have been gained by means of computer simulations 
@, 0, Iq-tlOJ. However, notwithstanding the significance of translocation phenomena both as possible technological 
application and from the standpoint of basic research, the understanding of the polymer translocation through a 
narrow pore is still elusive and in some respects controversial |ll| . 

In most theoretical descriptions (except Refs. [Tljj , where a full Cartesian-space treatment has been suggested) the 
translocation process is assumed to be captured by a single variable, the translocation coordinate s(t), which measures 

J" ' ■ the number of translocated segments at time t. Generally, s(t) would depend on the external potential, if present, 
and be influenced by random noise. At present it is well established that the translocation coordinate s(t) follows an 

^C) • anomalous, sub-diffusive law, although the physical background as well as the equation which governs the probability 
distribution function (PDF) W(s,t) that s segments have passed through the pore at time t are still controversial. 

It was suggested recently [5[ that the translocation dynamics is governed by a fractional diffusion equation (FDE) 
with an (anomalous) diffusion exponent a < 1. Thus a description of the translocation process in terms of W(s, t) was 
derived [6| which appeared to agree favorably with Monte Carlo simulation results. This approach was generalized 
[7| to the case of biased translocation, driven by external potential, which was described by an appropriate fractional 
Fokker-Planck equation (FFPE). As for the non-driven translocation, however, the FFPE approach yields a non- 
Gaussian distribution W(s,t) whose width A(i), diverges sub-diffusionally with the elapsed time, A(i) oc t a . While 
this latter property may readily be verified in numeric experiments, the unambiguous determination of the precise 
functional shape of W(s, t) is rather difficult due to progressively deteriorating statistics of the distribution at late 
times. 

Meanwhile, several recent publications |13 - [l4| ■ devoted to the unbiased translocation dynamics of a Gaussian ID 



I 

a 
a 
o 
_o 

> 

oo 
en 

r^ 

o 
o 

> 

X 



chain [12[ as well as to that of a 2D self-avoiding chain [13|, have validated the sub-diffusive behavior of (s 2 (t)} with 



an exponent a ~ 0.8. Nonetheless, these new findings cast serious doubts as to whether the FFPE provides indeed 
an adequate description of nondriven translocation dynamics: 



• it was found [12Hl4j by means of computer simulations that the probability distribution W(s, t) of the translo- 
cation coordinate s, stays Gaussian for different time moments not exceeding the mean translocation time (r). 
At larger times this distribution attains a more complex form. 

• the long time tail of the first-passage-time (FPT) distribution was found to be exponentially decreasing thus 
challenging the power-law behavior, suggested earlier [6] within the framework of the FFPE. 



These findings question the validity of the FDE approach in the description of translocation dynamics. Indeed, one 
cannot reconcile the aforementioned controversial features of the translocation dynamics within a FFPE description. 

In the present work we revisit the problem and demonstrate that a proper Fokker- Planck equation of motion (FPEM) 
which governs the PDF W(s, t) and faithfully reproduces all recently found properties of translocation dynamics may 
be rationalized and solved within the framework of fractional Brownian motion (fBm) [15]. Namely, we treat the 
definition of the translocation coordinate velocity, v(t) — ds(t)/dt, as a Langevin equation, where the velocity auto- 
correlation function, G(t) = (w(i)u(O)), is considered as dynamic input without making any a priori assumption about 
the dynamics that underlies the translocation process. In principle, G{t) may be taken from the computer simulation. 
Such an approach finally leads to a FPEM (with a time-dependent diffusion coefficient) for W(s, t) in unbiased and 
biased cases. We make an extensive Brownian dynamics simulation study (using GROMACS simulation package) in 
order to check and justify our analytical results. 

In Section |TI] we derive a governing equation for PDF W(s,t), calculate its first moments, and present an exact 
analytic solution for the first passage time distribution of translocation times as well as for the so called survival 
probability in terms of monomer consecutive number s and elapsed time t. A comparison with our Brownian Dynamics 
simulation data is performed in Section IIIII Then our main conclusions are given as a brief summary in Section IIVI 



II. CHAIN TRANSLOCATION DYNAMICS 
A. Fokker-Planck equation with a time-dependent drift and diffusion coefficient 

One may derive the FPEM for the distribution W(s,t), starting from the Langevin equation, 



dt 



s{t) = v{t), 



(1) 



where, by assumption, the translocation velocity v(t) = ds(t)/dt follows Gaussian statistics (i.e., is Maxwell-Boltzmann 
distributed). If a (generally time dependent) external driving force /(£) is present, the mean velocity (v(t)) = /(i)/£o 
where £o denotes the friction coefficient. Note that the velocity v(t) in eq. ([1]) reflects the change in the s coordinate 
per unit time. By measuring s(t) as the contour length of the chain on the trans-side of the membrane, a reasonable 
estimate for the translocation velocity is obtained. In fact, v(t) can be approximated by the velocity of the bead inside 
the pore in the direction perpendicular to the wall, which we denote by v z (t) (see Fig. [I]). Thus we tacitly assume 
that translocation velocity may be faithfully characterized by bead velocity in Cartesian space. The correctness of 
this approximation will be addressed in Section IIIII 




TRANS 



FIG. 1: (Color online) Chain translocation through a nanopore with an instantaneous translocation coordinate s(t). The 
translocation velocity is approximated by the velocity of the bead which resides inside the pore. The separating membrane is 
taken as sufficiently thin so that no more than one bead may be in the pore at a time. The driving force caused by the chemical 
potential gradient within the pore, / = (/ii — ^2)/T, acts on the same bead, /ii and /i2 denote values of the chemical potential 
on both sides of the separating membrane. 

The variance of s(t) can easily be calculated under the assumption that the velocity auto-correlation function 

(v(ti)v (£2)), defined as (v(ii)u(£ 2 )) = limT^oo(l/T) J v(t\ + t)v{t 2 + t)dt depends on the difference |£i — t 2 \ only 
(assuming time translational invariance). This leads to the relationship [161 ] 



({s- So } 2 (t))=2 / (t-t'){v(t')v(0))dt'. 
Jo 



(2) 



where so = s(t = 0) is the initial value of the translocation coordinate. In the case of biased translocation, the driving 
force is usually due to a chemical potential gradient A/i, which is typically generated by applying a voltage difference 
across the membrane, i.e., / = A fi/T, where A fj, = /j,\ — fi2- Here / is imposed on the monomer which is currently 
in the pore (cf. Fig. Q}. As pointed out above, the velocity v(i) is assumed to be a Gaussian random variable with 
the first two cumulants given as 

(v(t)) = f]p-, (3) 

and 

G(h,t 2 ) = ([v(h) - (v(h))Mt 2 ) - (v(h)))). (4) 

We next consider the PDF W(s,t), which is defined as 

( 
W(s,t) = (s(a - so - J v{t')dt')), (5) 

o 
and satisfies the initial condition W(s, 0) = S(s — sq). It is straightforward to show that 
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° W (s, t) = -£- [(v(t)) (6(8 - so - / v(tf)dt'))] - i- / drG(t, t) (^- S(s - s - v{t')dt')), (6) 



ds J ' \5v(t) 





where we have used Novikov's theorem [T3, [l^l with the two first cumulants of v(t) given by eqs. <|3j) and (J4]). In the 

t t 

last term in eq. ((6|) we have also used ((8/8v(t)) 8(s — sq — J v(t')dt')) — —(d/ds)(8(s — sq — J v(t')dt')) as well as 



Jo S(t - t')dt' = 1 for t < t. Eventually the FPEM takes on the form: 

§-w(s lt) = -m°w(s lt ) +m -f 

ot £0 OS OS 



W(s,t) = -^ — W{s,t) +V{t)—r W{s,t). (7) 



(■I 



with a time-dependent diffusion coefficient T>(i) = L G(t,r)dT. This equation of motion is different from the FFPE 
19] which was used in our previous investigation [3 0|- In contrast to FFPE, it is Markovian but contains a time- 
dependent drift term and a time-dependent diffusion coefficient. Moreover, the form of eq .([7]) does not depend on 
the underlying dynamics (underdamped or overdamped regimes!) as far as the driving force f(t) and the velocity 
correlator G(ti,t2) are not specified. We should emphasize that Eq. ([7]) has been by assuming that the effective 
friction coefficient £n r emains constant with time. This is supported by computer simulation data, cf. Fig. 6a in 
Bhattachrya et al. [201 ]. I was shown [2fj] that for a driven translocation the mean center-of-mass velocity quickly 
saturates to a constant value after a short initial increase. The transient regime itself lasts less than 10% of the total 
translocation time for a short chain (N — 128) and quickly vanishes for longer chains. Recently the translocation 
dynamics has been treated [2l|, [22| on the basis of a FPEM with a time dependent drift (even for a constant driving 
force) and diffusion coefficient. 

In the simplest case of an unbiased process (i.e., / = 0), the system is in a steady state and T>(t) = L G(t — r)dr = 

J Q G(r)dr. Then with the boundary conditions at infinity eq.© may be solved by defining a new time variable 



t 
f t = Jdt'G(t% (8) 





t t' t 

i.e., 9(t) = Jdt' J dt" G{t") = j(t - t')G(t')dt'. A comparison with eq.© then yields 
00 

6{t)= l -{(s-so) 2 )<xt a , (9) 

where we have assumed that the sub-diffusive behavior of ([s — s ] 2 (t)) is described by an exponent a < 1. Thus eq. 
© takes on the form 

^W(s,e) = ^W(s,9). (10) 



The solution of eq. <jTU]) reads W(s, 9) = (1/2tt9) 1 ^ 2 exp [— (s - s ) 2 /4#] . Then, after taking into account eq.©, one 
finds 



W(s,t) 



1 



VZnDt 



exp 



(g ~ go)' 
ADt a 



(11) 



where O is a constant. This result for W(s, t) has been obtained very recently by Panja |23l . 124} in a more general 
context dedicated to anomalous polymer dynamics. 

Evidently, the distribution Wis, t) is indeed Gaussian albeit with a width proportional to the second moment that 
grows sub-diffusively with elapsed time. This result reproduces the Monte-Carlo simulation findings, reported recently 
by Kantor and Kardar [l2|, L!j|. It should be noted that such a process, described by a Gaussian distribution with 
anomalous width is sometimes referred to as fractional Brownian motion |15l 125} and has been discussed in the context 
of DNA sequence statistics [26| . 

B. Boundary conditions 

Eq. (|10j) should be considered along with the appropriate boundary conditions. In the special case when the 
translocation starts from the middle monomer (i.e., sq — N/2) and stops as soon as any of the two chain ends passes 
the pore, both boundary conditions are adsorbing, i.e., W(s — 0,9) — 0, W(s = N,9) — 0. The corresponding initial 
condition reads W(s, 6 = 0) = S(s - s ). Recall that 9{t = 0) = 0. 

The full solution can be represented as a sum over eigenfunctions (fi n (s), i.e., W(s,9) — Y^=oT n {9)ip n (s) (see, 
e -g-i [27} ) where f n (s) obey the equations [d 2 / ds 2 )ip n (s) + \ n ip n (s) = 0. The temporal part is governed by the 
equation (d/d9)T n {9) + \ n T n {9) = with eigenvalues A„. The proper eigenfunctions which satisfy both boundary 
conditions, are (p n (s) = A n sin(n7r/iV) (the factor A n can be fixed by the initial condition). The eigenvalues are then 
A n = {im/N) 2 . The completeness relation for the eigenfunctions (see Sec. 5.4 in [23] ) has in this case the form 

2 P®_ / T17TS \ / TITTSn \ 

5(s — so) = — V sin ( ) sin ( ) so that the full solution becomes thus 

^V n =o \ N J V N ) 



in*. (h = — Y, sin i^jf) sin (^r) ex p 



N 



n=0 



/nir\- 

v"/vJ 



(12) 



Note that eq. (|9]) defines a sub-diffusional law 9(t) = D t a , which governs the solution, eq. (fl"2"]) . Eventually, if the 
starting point so = N/2 (as is the case in our MD-simulation) , only the odd terms in the series, eq. (|12[) . survive, i.e., 
n = 2m + 1 and the final result reads: 



W(s,t) = - £(-!)" 



sm 



(2m + l)7rs 



N 



exp 



(2to+1) 2 7t 5 

N~ 2 



Dt a 



(13) 



The first two moments generated by this probability distribution contain a lot of information, which permits us to 
compare our theoretical results with those from literature and with the numerical results of Section IIIII 



C. First and second moments of W(s,t) 



Taking into account that J Q sin[(2m + l)7rs/iV]ds = 2N/(2m+l)ir and J Q ssm[(2m + l)irs/N]ds = N 2 /(2m+ l)n, 

N N 

one obtains for the first moment (s) = J sW(s,t)ds/ J W(s,t)ds = sq = N/2. The second moment is centered at so, 

o o 

N N 

hence, A = (s 2 ) — (s) 2 = j [s — so] 2 W(s,t)ds/ J W(s,t)ds. This leads to the following expression for the variance, 

o " o 



N 2 



oo C 1 \m 

m= o (2TO+1) 3 


N 2 


a ^ (-l) m 

7T > , -, exp 

rto (2m + 1) 


(2m+l) 2 TT 2 1 
v J Dt a 

N 2 



(14) 



At t — >• oo the second moment A reaches a plateau, i.e., A = (7V 2 /4)(1 — 8/ir 2 ). It can be readily shown (taking 
into account that Em=o(-l)"7( 2m + 1) = V 4 and Em=o(- 1 ) m /( 2m + I) 3 = ^ 3 /32) that A(£ = 0) = 0. The time 
dependence of the second moment A, given by eq. (TT4")) . will be discussed further in the Section ITO1 



D. First passage time distribution 



The first passage time distribution (FPTD) function is defined as [27J and describes the probability distribution of 
the observed translocation times: 



JV 

Q{t) = -- J W(s,t)ds. 



(15) 



Taking into account eq. (I13p . this definition yields 



Q(t) = 



N 2 



Y^ (-l) m (2m+l) exp 



m=0 



2^2 



(2m + 1) 2 tt 

w 



■Dt a 



Q(t) 



AiraDt 



n* ™n-w Dta 



(16) 



The long time behavior of Q is determined by the smallest eigenvalue in the series eq. (flB]) . Therefore, at t — > oo one 
has 



(17) 



i.e., the FPTD follows a stretched-exponential law at late times. We will try to check this prediction in our MD- 
simulation study in Section [IIII 



E. Asymptotic behavior near the adsorbing boundary 



An interesting aspects of translocation dynamics in a system with two adsorbing boundaries has been considered 
recently [12l4l4l |. suggesting that for a sufficiently long time interval the normalized distribution 



Ps (s,t) 



W(s,8 ,t) 

f^W{s,s ,t)ds' 



(18) 



reaches a stable (time- independent) shape which differs from the simple sin- function and at s — > (i.e., close to the 
adsorbing boundary) the functio n p ? n (s) ~ s^ where 4> > 1 (for a sin-function it would have been <f> = 1). 

Based on numerical results [ijg, [Da] > one nas tried [2^| recently to link translocation dynamics to self-affine processes 
undergoing anomalous diffusion in bounded domains within the context of fBm. The argumentation [281 ] relies on two 
crucial points: (i) The PDF W(s,s ,t) has a self-affine form W(s,s ,t) = (l/t a / 2 )F(s/t a / 2 ,s /t a / 2 ), and (ii) the so 
called survival probability S(sq, t) = L W(s, so,t)ds has a long-time scaling behavior S(sq, t) ~ t~ e with an exponent 
= 1 — a/2. Simple scaling arguments lead Zoia et. al. 28( to the conclusion that the normalized probability 
p So (s,t) = {\/t a / 2 )py {y), where y = s/t a / 2 , j/q = s o/t a ^ 2 and the scaling function p yo (y) ~ y^, with the exponent 
(j) = 2/a-l> 1. 

Note that our solution, eq. (fTTj) , disagrees with the scaling proposed in [28j. Since the eigenvalue spectrum is 
discrete (as it should be for a finite interval < s < N), only the smallest eigenvalue dominates in the long time limit 
so that for the PDF W(s, So, t) and the survival probability S(so,i) we get 



W(s,s ,t) 
S(s ,t) 



_j sul (_j sln (_j exp 



(-) Dt" 



%**m->-{v)>*- 



Therefore, 



Ps (s,t) 



W{s,s ,t) 
S(s ,t) 



7T /7TS 

sin — 

27V \JV 



(19) 



(20) 



Evidently, the self-affine scaling [28], conjectured for W(s, So,t) and S(so,t), does not hold and the "stable shape" 
at t — > oo is a simple sin- function. This behavior is in accordance with the definition of fBm [Tjj. The stretched 
exponential behavior for the survival probability, which corresponds to eq. I J19I ), has been recently discussed in ref. 



III. SIMULATION RESULTS 

To study the translocation process numerically, we performed Brownian Dynamics (BD) simulations for two different 
chain lengths : N = 51 and N = 101. The polymer was modeled using a coarse-grained description in which the 
adjacent monomers are connected by finitely extensible nonlinear elastic (FENE) springs, corresponding to a pair 
potential 



U 



FENEVij 



(r, 



k/i 



In 1 



Rl 



(21) 



where r-y is the bond length between two beads and Rq = 1.5 denotes its maximal extension. All beads experience 
excluded volume interactions which are modeled by a Lennard- Jones potential Ulj, defined by 



EMr«)=4e 



12 



(22) 



where we use a cut-off r c = 2~ 1 ' 6 er, implying that Ulj = for r^ > r c . The parameter values were taken as 
e = 1.0,cr = 1.0, k = 30.0, and were kept fixed during the simulations. The friction parameter £ was taken as 
£ = 100-y/me/cr, and the temperature T — 1.2e/ks, this implies a monomer diffusion coefficient Dq = fc,eT/£ = 1.2e/£. 

In order to simulate a translocation event, we create a separating membrane, which consists of monoatomic layer, 
with a hole (the center monomer removed). The membrane is placed inside a box of size (54x54x54) with periodic 
boundary conditions and lies in the z = 0-plane. All atoms in the membrane are frozen and interact with the 
translocating chain via the Lennard- Jones potential. Before the translocation has started, we put the middle bead of 
the chain symmetrically inside the pore, and equilibrate the configuration, while keeping the middle bead fixed. To 
verify that the chain has completely relaxed to equilibrium, we wait until the average radius of gyration no longer 
changes with time. We have checked a posteriori that the number of translocation events to the trans and cis side 
are evenly distributed. 

The translocation process of the polymer chain is studied by Brownian dynamics simulations in a fictitious solvent. 
The translocation velocity v(t) is sometimes approximated by v z (t), the z— component of the velocity of the bead 
inside the pore. From the recorded data of the translocation coordinate and the velocity v z (t) we make histograms for 
the velocity- and the s— coordinate distributions; see also Figs. [3] and |3^. It is convenient to define the translocation 
coordinate s(t) as a continuous variable so that one may easily calculate the translocation velocity v(i) = ds(t)/dt. In 
Fig. [2] we show a simple way to do this. The origin of the coordinate system is located in the pore. The z-coordinates 
of the n-th and n + 1-th beads are marked as z n and z n +i respectively so that z n > and Z n +x < 0. The continuous 
translocation coordinate which interpolates between two consecutive integer values n and n + 1 can be defined as 




n+1 



, if z n ^ and z n+1 ^ 0, 
,if z n = 0. 
,if z n+ i = 0. 



(23) 




FIG. 2: (Color online) The individual chain segments are enumerated starting from the right-most terminal bead. The origin 
of the coordinate system is placed in the pore. The z-coordinates of n and n + 1 beads are marked as z n and z n +i respectively. 



We plot the the velocity distribution function P(v z ) (see Fig.|3]and the translocation coordinate distribution W(s, t) 
(Fig. |1 a) for several times t = 1.0 x 10 4 , t = 5.0 x 10 4 , t = 1.0 X 10 5 , t = 1.5 X 10 5 , and t = 5.0 X 10 5 time steps. 



All distributio ns w ere obtained by averaging over at least 5000 runs. Note that 1 time step corresponds to 0.005 
ps = 0.005a^m/e. The velocity distribution is clearly seen from Fig. [3] to be a Gaussian, centered around v z = 0. 
We also plotted P{v z ) for t — 50000 and t = 500000 time steps so as to demonstrate that the distribution is indeed 
time-independent. Moreover, we verified that the distribution is well described by a Maxwell-Boltzmann distribution 
for the simulation temperature T = l.2e/ks- The PDF for the s— coordinate, W(s,t), is centered around s = 51, 
which is the middle bead of the chain of length N — 101; the distribution broadens symmetrically. From Figs. [3] and 
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FIG. 3: (Color online) The velocity distribution function obtained from numerical simulations at T — 1.2e/fes, averaged over 
9000 runs for two different times (t=250 ps and t=2500 ps). The simulations all collapse on a single Maxwell-Boltzmann 
distribution curve plotted in red. 



0k, it can easily be seen that the distributions for both the velocity and the translocation coordinate s, obtained from 
our Molecular Dynamics simulation are indeed Gaussian. In Fig. 2)3 we show change of the variance A(i) (as defined 
by eq. (Tl4l) ) of the Gaussian distribution curves W(s, t) as a function of time. 
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FIG. 4: (Color online) (a) The PDF of the translocation coordinate W(s,t) at 5 different time moments. The 5 Gaussian fits 
correspond to the PDF of s at time steps 10 4 till 10 6 and have a variance A(t) whose value is given in the legend and which 
increases with the number of time steps, (b) The variance A(t) of the distribution W(s, t) as a function of number of time 
steps. It is clear that the behavior for short times is sub-diffusive and scales as t 0,55 . For larger times t > 10 5 time steps, the 
behavior changes and is only slightly sub-diffusive and A(i) ~ i ' 91 



Evidently, A(t) behaves sub-diffusively as ex i 055 for shorter times, and later switches to a r°- 91 -law for longer 
times. This corresponds to the simple scaling consideration [30] which claims that the variance of the translocation 
coordinate A(i) = (s 2 ) — (s) 2 in the long time limit should go as A(t) oc i Q , where a = 2/(2i/+ 1) sa 0.92 for 3D-case. 



8 



This power is obtained by the requirement that the translocation time Tt r should scale like the maximal relaxation 
time of the polymer chain, i.e., T tr ex N 2u+1 . In addition, by the time when the translocation is completed, the 
following relationship holds: A(rt r ) oc iV 2 . 

We emphasize that the change in the exponent a from 0.55 in the range [10 4 , 10 5 ] time steps to 0.91 from [10 5 , 10 6 ] 
time steps, is in agreement with recently reported results by Amatai et. al. [14j . who considered the movement of a 
tagged monomer in a Rouse chain. This somewhat surprising result indicates that the study of the movement of a 
tagged Rouse monomer can indeed provide essential information about the translocation process. 

Next we consider the velocity correlation functions. The normalized auto-correlation function is depicted in Fig.[5ji. 
From this figure one can infer that the correlation time is not zero albeit rather small. One clearly sees anti- 
correlation in the interval t 6 (0,5), after which the correlations vanish on the average. Note that such unusual 
'negative' correlations at short time lags have been observed recently in experimental studies of the sub-diffusive 
motion of bacterial chromosomal loci through a viscoelastic cytoplasm 3l| . To illustrate the invariance with respect 



to shifts in time, a time shifted correlation function is also displayed. Since the curves collapse onto each other, this 
demonstrates that the correlation function is indeed time invariant as required. In Fig. [5}d we give the consistency 
check: the translocation coordinate mean-square displacement A(t) as derived from the MD-simulation is compared 
to the time integral given by eq.©. One can see very good agreement between these two ways of A(t) calculation, 
which also suggests that the velocity correlator (cf. Fig. [5£i) has been found with reasonable accuracy. 

We have also relaxed the time translational invariance assumption and represented A(i) as double integral 

J dt' J Q dt" (v(t')v(t")) . The results of these calculations basically superimpose onto each other which indicates 
again that the time translational invariance holds indeed. From Fig. \5\ it can also be seen that the z-component of 
the monomer velocity inside the pore provides a reasonable measure for the translocation velocity although there 
definitely are minor deviations between this velocity and the translocation velocity. In particular, the variance of the 
tagged particle velocity has a somewhat smaller power law exponent than the real translocation velocity. However, 
combining Fig. [5Jd with Fig. |4Jd shows that the exponent a changes with time from 0.91 on the very short time scales 
[1, 200] to about 0.55 for times on the order of 10 4 . This is again in agreement with the prediction of [14[ for a tagged 
Rouse monomer. 
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FIG. 5: (Color online) (a) The normalized velocity correlator against elapsed time t. (b) Comparison of the A(t) data 
(circles) with the result of integration of the velocity correlator according to eq.© (blue solid line). The double time integral 

L dt' f dt" (v(t')v(t")), which is valid even when time translational invariance does not hold, is given by green + symbols. 



Figure |5] demonstrates the A(i)-behavior for a larger time interval where it is seen to reach a plateau as soon as the 
polymer chain has completed the translocation process. Here we compare the variance A(i) from numerical simulation 



with the analytical expression eq. (fT4|) where we took a — 0.92. It can be clearly seen that the results are in good 
agreement. 
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FIG. 6: (Color online) (A(t)) for long times reaches a plateau. The plateau height as well as the slope of the theoretical curve 
preceding the plateau are in a good agreement with simulation results. 



The first passage time distribution Q(t) is interesting too. In a previous paper [6j, we employed the FFPE to 
describe the distribution Q(t). The tail of the distribution decreased as a power law in that case. However, for the 
theory that we develop here, a stretched exponential decay is expected. For a N — 51 chain, we have performed 
simulations to compare Q(t) with the theoretical prediction. For an average over 400 runs, we obtained the result 
which is shown in Fig. [7J The theoretical distribution is found to be in reasonable agreement with the simulation 
result. We remark that experimental and numerical verification of the tail of the translocation time distribution is 
extremely difficult, as the events constituting the distribution's tail are very sparse. 
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FIG. 7: (Color online) The FPT distribution Q(r). Theory and simulations compared for a N=51 chain. 



IV. CONCLUSION 



In the present work we treat the translocation dynamics within the framework of fractional Brownian motion 
(fBm) whereby the definition of the translocation coordinate velocity, ds(t)/dt = v(t), is considered as a Langevin 
equation where v(t) can be viewed as a Gaussian process with given two first cumulants: (v(t)) = f(t)/^o and 
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G(t) — (v(t)v(0)) — (v(t)) 2 . The first cumulant implies that the friction £o is constant, in agreement with data 
from computer experiments. This does not rule out the possibility that there are cases when friction might be 
time-dependent [ljj. In such cases a modified version of the governing equation (|7|) would be warranted. 

Based on this Langevin equation, we derive (without any further assumption) the appropriate Fokkcr-Planck equa- 
tion of motion for the distribution function W(s, t) with a time-dependent drift and diffusion coefficient. The obtained 
solution for W(s, t) demonstrates two characteristic features of fBm which agree favorably with the recent findings 
[12M14J ] as regards nondriven translocation dynamics: 

• a Gaussian distribution for the translocation coordinate s during the translocation process 

• a sub-diffusive behavior for the variance of the translocation coordinate, A(£) = (s 2 (£)) — (s(t)} 2 . 

Moreover, the survival probability S(t, Sq) in the long time limit has a stretched-exponential form in agreement with 
recent findings [29| (albeit in contrast to a popular opinion about its scaling behavior [231). O ne should note, however, 
that the power (0.92) of the stretched exponential function is very close to unity so that a distinction from simple 
exponential decay of the survival probability at the long time limit would be hardly detectable. 

The velocity correlator G(t) for the case of unbiased translocation has been computed by means of Brownian 
dynamics simulationi and found to expose unusual negative correlations in time. A consistency check, implemented 
by means of eq.©, shows that the calculated G(t) is meaningful indeed. Thus, it appears that the presence of negative 
velocity correlations reveals a feature which is specific for the process of translocation in particular and for anomalous 
diffusion in general [3l| . 

We have calculated the variance of the translocation coordinate A(£) and shown that A(£) follows a power law, 
with different exponents a, depending on the timescales one is interested in. In particular, for times on the order 
of the translocation time, A ~ £ 091 in 3D. Furthermore, we examined the translocation time distribution function 
Q(t), derived from simulation results in a broad time interval and compared these to our theoretical predictions. 
Unfortunately, attaining good statistics for very late times still remains a difficult task for present days simulation 
studies so more work will be needed in order to test the agreement with theoretical predictions unambiguously. 

Eventually, we have shown in Section |TT] that our approach can be used for the biased translocation as well. This 
investigation will be reported on in a separate presentation. 
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